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ABSTRACT 

We analyze a sample of 32 galaxies for which a dynamical estimate of the mass of 
the hot stellar component, Mbuige; is available. For each of these galaxies, we calculate 
the mass of the central black hole, M,, using the tight empirical correlation between M, 
and the bulge stellar velocity dispersion. The frequency function A^[log(M,/Mbuige)] is 
^ \ reasonably well described as a Gaussian with (log(M,/Mbuige)) ~ —2.90 and standard 

\Q ' deviation ~ 0.45; the implied mean ratio of black hole to bulge mass is a factor ~ 5 

t^ ■ smaller than generally quoted in the literature. We comment on marginal evidence for 

Qv^ ■ a lower, average black-hole mass fraction in more massive galaxies, which should be 

^^ , investigated using larger samples. The total mass density in BHs in the local Universe 

f~^ ! is estimated to be ~ 5 x 10^ M© Mpc~^, consistent with that inferred from high 

redshift {z ~ 2) AGNs. 

o 



1. Introduction 

With an ever-increasing number of secure detections, supermassive black holes (BHs) have 
evolved, in a ten-year span, from exotic curiosities to fundamental components of galaxies. It is 
;h ' now generally accepted that the formation and evolution of galaxies and supermassive BHs are 

tightly intertwined, from the early phases of proto-galactic formation (Silk & Rees 1998; Sellwood 
&; Moore 1999), through hierarchical build-up in CDM-like cosmogonies (Efstathiou & Rees 1988; 
Haehnelt & Rees 1993; Haehnelt, Natarajan &: Rees 1998; Haiman &; Loeb 1998), to recent galaxy 
mergers (Merritt 2000). Studying the demographics of the local BH population might have a 
significant impact on models of galaxy evolution (e.g. Salucci et al. 1999; Cattaneo, Haehnelt & 
Rees 1999; Kauffmann & Haehnelt 2000). 

Magorrian et al. (1998) presented the first, and to date only, demographic study of nuclear 
BHs. Ground based kinematic data for 32 galaxies were combined with HST photometry 
to constrain dynamical models - based on the Jeans equation - under the assumptions of 
axisymmetry, velocity isotropy in the meridional plane and a spatially-constant mass-to-light ratio 
for the stars. The mass of a putative nuclear BH was introduced as a free parameter, in addition 
to the stellar mass-to-light ratio and the galaxy inclination angle. In most of the galaxies, the 
addition of a central point mass improved the fit to the observed kinematics. Magorrian et al. 



concluded that most, or all, galaxies contain central BHs with an average ratio of BH mass to 

ilgc 



spheroid mass of M,/Mbuigc ~ 10 ^. 



The Magorrian et al. study remains unique for targeting a large sample of galaxies, and for its 
coherent and homogeneous treatment of the data. However, while the Magorrian et al. estimates 
of the bulge mass-to-light ratios are likely to be robust, a number of authors have noted that the 
inferred BH masses might be systematically too large. Van der Marel (1997) showed that the BH 
masses derived from well-resolved central kinematical data are a factor 5 smaller than produced 
by the Magorrian et al. analysis; he suggested that the neglect of velocity anisotropy might 
have led to overestimates of the BH masses. Wandel (1999) compared BH masses derived from 
reverberation mapping studies of active galaxies with the Magorrian et al. estimates and found a 
discrepancy of a factor of ~ 20 in the BH-to-bulge mass ratio at a fixed luminosity. He noted the 
difficulty of resolving low-mass BHs in distant galaxies and suggested a distance-dependent bias 
in the estimates. 

An independent argument along the same lines was presented by Ferrarese &: Merritt (2000; 
FMOO; Paper I). Using the tight empirical correlation between M, and a^ the velocity dispersion 
of the stellar bulge, for the 12 galaxies with the best-determined BH masses, FMOO showed that 
the Magorrian et al. masses fall systematically above the M, — a relation, some by as much as two 
orders of magnitude. 

At the present time, the M, — a relation is probably our best guide to BH demographics. 
Ferrarese &: Merritt (2000) found that the relation has a scatter no larger than that expected on 
the basis of measurement errors alone. The relation is apparently so tight that it surpasses in 
predictive accuracy what can be achieved from detailed dynamical modeling of stellar kinematical 
data in most galaxies. By combining the bulge stellar masses derived by Magorrian et al. with 
BH masses inferred from the M, — a relation, we are in a position to compute the most robust 
estimate to date of the BH mass distribution in nearby galaxies. 



2. Data 

Table 1 gives the relevant physical parameters for the 32 galaxies in the Magorrian et al. 
sample. All galaxies, with the exception of M31, are early type. In what follows, we refer to the 
hot stellar component in these galaxies as the "bulge;" this is in fact the case for M31, although for 
the other objects the "bulge" is the entire galaxy. Distances were re-derived as in Paper I; values 
for the bulge F-band luminosity (Lbuige), bulge mass (Mbuige) and BH mass (Mfit) are the same 
as in Magorrian et al. except for the (mostly small) corrections resulting from the new distances. 

Central velocity dispersions Oc were taken from the literature and corrected to a common 
aperture size of 1/8 of the effective radius, as in Paper I. We then computed BH masses, M,, using 



the M, — (Tc relation in the form given by Merritt & Ferrarese (2000; Paper II): 

M, = 1.30 X 10^ Mo(o-c/200 km s"^)^-^^. (1) 

This expression was derived by fitting to the combined galaxy samples of Ferrarese & Merritt 
(2000) (12 galaxies) and Gebhardt et al. (2000a) (15 additional galaxies), plus 7 active galaxies 
for which both a^ and M, are available, the latter from reverberation mapping (Gebhardt et al. 
2000b). Some debate exists over the exact value of the slope in equation (1) (Paper I, II; Gebhardt 
et al. 2000a). We explore below how changing the assumed slope affects our conclusions. 

The correlations between M, and Lbuigo and between M, and Mbuigo are shown in Figure 1. 
There is a rough proportionality of both Lbulgc and Mbuigc with M,, though the vertical scatter in 
both relations is much larger than in the M, — a relation (Paper I). 

We defined the two mass ratios: 

Xfit = Mfit/Mbulge, 

X = M./Mbulge (2) 

based respectively on the BH mass estimates from Magorrian et al. and from the M, — a relation. 
Values of logxgt and logx are given in Table 1. BH masses derived from the M, — a relation 
yield the mean values {x) = 2.50 x 10^'^ and (logx) = —2.90. These are substantially smaller 
than the mean values computed from the Magorrian et al. BH masses: (xfit) = 1.68 x 10^^ and 
(logXfit) = —2.20. We note that one galaxy, NGC 4486b, has logxgt = —0.54, making it an 
extreme outlier in the Magorrian et al. distribution. Removing this single galaxy from the sample 
gives (xfit) = 7.2 x 10~^ while leaving (logXfit) essentially unchanged. 

Figure 2 reveals a clear trend of Mfn/M, with the apparent radius of influence of the 
central black hole, assuming the masses predicted by the M, — a relation are correct. A natural 
interpretation is that there is a resolution-dependent bias in the Magorrian et al. modeling (e.g. 
van der Marel 1997; Wandel 1999): the radius of influence of all of the Magorrian et al. galaxies is 
smaller than 1 arcsec, too small to have been clearly resolved from the ground. 



3. Analysis 

We seek an estimate of the frequency function N{y) = N(\ogx). Following Merritt (1997), we 
define this estimate as N{y), the function that maximizes the penalized log likelihood 



log Cp = J2 log(^ ° E), - XP{N) (3) 

of the data yi,i = 1, ..., n, subject to the constraints 

N{y)dy = 1, N{y) > 0. (4) 



Here N o E is the "observable" function, i.e. the convolution of the true N with the error 
distribution of y. This error distribution is not well known; we assume that it is a Gaussian with 
some dispersion A^^. Failing to account for measurement errors in y would lead to a spuriously 
broad N{y). 

The natural penalty function to use is Silverman's (1982): 

P{N)= / {d/dyf log N{y) dy. (5) 

This function assigns zero penalty to any N{y) that is Gaussian. In the limit of large A, the estimate 
N is driven toward the Gaussian function that is most consistent, in a maximum-likelihood sense, 
with the data; smaller values of A return nonparametric estimates of N{y). We made no attempt 
to calculate the "optimal" value of the smoothing parameter given the small size of the data set. 

The results of the optimization are shown in Figure 3, assuming Aj^ = 0.15. N{y) is nicely 
symmetric and reasonably well described as a Gaussian, although with a narrower-than-Gaussian 
central peak. The best-fit Gaussian has its mean at y = logx = —2.93 and a standard deviation 
of 0.45. 

By contrast, the Magorrian et al. masses define a more flat-topped distribution with one 
extreme outlier, NGC 4486b, at logXfit = —0.54. The Gaussian fit to the Magorrian et al. mass 
distribution has its mode at —2.25 and a standard deviation of 0.52. 

The two galaxies with the largest BH mass ratios, NGC 4486b and NGC 4660, are both 
low-mass ellipticals. The smallest mass ratio, logx = —3.95, is seen in a very massive galaxy, 
NGC 4874. It is therefore interesting to check whether low- and high-mass galaxies have different 
characteristic distributions of logx. This hypothesis is tested in Figure 4a, which shows N{y) 
computed separately for the 16 galaxies from Table 1 with the lowest and highest values of 
-^buige- There is in fact a slight difference between the two distributions: the high- mass galaxies 
have (logx) = —3.10 and a\ogx = 0.39, while the low-mass galaxies have (logx) = —2.71 and 
o'logx = 0.49. However the offset in (logx) is similar to the width of either distribution and may 
not be significant. We note that the massive galaxies define a narrower distribution. 

Our conclusions about N{y) might be substantially dependent on the assumed form of the 
M, — a relation, equation (1). The slope of that relation is fairly uncertain, a = 4.72 it 0.4; 
however the normalization at cjc ~ 200 km s~^appears to be more robust (Paper II). We therefore 
set 

M, = 1.30 X 10^ M0(o-c/2OO km s'^)" (6) 

and investigated the effects of varying a. Figure 4b shows that a larger slope implies a broader 
N{y) due to the stronger implied dependence of M, on ac- However the mean value of logx is 
almost unchanged. 



4. Discussion 

Our estimate of the mean BH-to-bulge mass ratio, (logx) ~ —2.90, falls squarely between 
the estimates of Magorrian et al. (1998) (~ —2.28), based on dynamical modeling of the same 
sample of galaxies used here; and of Wandel (1999) (^ —3.50), based on BH masses computed 
from reverberation mapping in a sample of 18 active galaxies. 

Bulge masses in the Wandel (1999) study were computed directly from bulge luminosities 
assuming a simple scaling law for the mass-to-light ratio, and not from dynamical modeling. 
There is reason to believe that these luminosities are systematically too large and therefore that 
the derived mass ratios M,/Mbuige are too low. Gebhardt et al. (2000b) and Merritt & Ferrarese 
(2000) found that the reverberation mapping BH masses in 7 galaxies were consistent with the 
M, — a relation even though they fall systematically below the M, — Lbuigc relation. A reasonable 
conclusion is that the true or derived luminosities of these active galaxies are systematically higher 
than those of normal galaxies with comparable velocity dispersions. A mean offset of a factor ~ 4 
in the bulge luminosities would suffice to bring the average mass ratio for active galaxies in line 
with the value inferred here. Gebhardt et al. (2000b) discuss a number of possible reasons why an 
error of this sort is likely in the AGN bulge luminosities. 

The discrepancy with the Magorrian et al. (1998) masses is perhaps unsurprising given past 
indications that these masses are systematically too large (van der Marel 1997; Ho 1999). The 
difference between (logx) and (logxgt) corresponds to a factor ~ 5 average error in the Magorrian 
et al. BH masses. One possible explanation is the neglect of anisotropy in the modeling (van der 
Marel 1997), but we emphasize that the errors in Mgt implied by Figure 2 are enormous, of order 
10 — 100, in many of the galaxies. If the BH masses predicted by the M, — a relation are correct, 
the kinematical data for these galaxies would not have contained any useful information about the 
mass of the BH (Figure 2). Any features in these data that were reproduced by adjusting Mgt 
must have had their origin in some systematic difference between the models and real galaxies, 
and not in the gravitational attraction of the BH. This conclusion, if correct, underscores the 
dangers of an "assembly-line" approach to galaxy modeling. 

We may crudely estimate the total mass density of BHs in the local universe by combining 
our result, M,/Mbuige ~ 1-3 x 10^^, with the mean mass density of spheroids, Pbuigc ~ 3.7 x 10^ 
M0Mpc^^ (Fukugita, Hogan & Peebles 1997, for Hq = 75 km s^^Mpc^^). This simple argument 
(first invoked by Haehnelt, Natarajan & Rees 1998) gives p,^L ~ 4.9 x 10^ MqMpc"^. Salucci et al. 
(1999) presented a more sophisticated treatment based on convolution of the spheroid luminosity 
function with iV(logx). They assumed a Gaussian distribution with (logx) = —2.60 and found 
p,^L --^ 1.7 X 10^ MqMpc^^. Correcting their value of (logx) to our value of —2.90 implies a factor 
~ 2 decrease in /9,,l, consistent with the result of our simpler calculation. 

The total mass density of BH at large redshifts can be estimated using an argument first 
suggested by Soltan (1982). Requiring the optical QSO luminosity function to be reproduced 
purely by accretion onto nuclear BHs, and assuming an accretion efficiency of 10%, leads to 



/9,,^ ~ 2 X 10^ MoMpc^^ (Chokshi & Turner 1992; Salucci et al. 1999). While independent of the 
cosmological model, this result is subject to uncertainties in the bolonietric corrections applied 
to the QSO magnitudes (e.g. Salucci et al. 1999), furthermore, concerns have been raised about 
the completeness of the QSO luminosity function (e.g. Goldschmidt & Miller 1998; Graham, 
Clowes & Campusano 1999). A similar argument, based on the hard X-ray background, gives 
P»,z ~ 3 — 4 X 10^ M0Mpc~'^ at z ~ 1.5 (Fabian h Iwasawa 1999; Salucci et al. 1999; Barger et al. 
2000). These numbers are consistent with our estimate of Pm,L- 

By contrast, p,^^ differs from the local BH mass density implied by the Magorrian et al. 
relation by over an order of magnitude, assuming a canonical 10% accretion efficiency onto the 
central black hole in high redshift AGNs. Haehnelt, Natarajan &: Rees (1998) and Barger et al. 
(2000) point out that if the remnants of the QSOs are to be identified with the BHs in present-day 
galaxies, the Magorrian et al. mass distribution requires either that a large fraction of BHs reside 
within high redshift sources that are too obscured (both in the optical and the X-rays) to be 
observed, or else that a significant amount of accretion (with low radiative efficiency) proceeds to 
the present epoch. The need for these alternative explanations is largely removed when the more 
robust estimate of p,^i^ presented in this paper is adopted. 
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Table 1. 



Galaxy 


Distance 


-^bulgc 


-'Wbulge 


l^c 


Mflt 


log ajflt 


M. 


logx 


N221 


0.8±0.1 


0.0373 


0.00807 


76±10 


0.0228 


-2.53 


0.0135 


-2.78 


N224 


0.8±0.0 


0.724 


0.350 


112±15 


0.598 


-2.75 


0.0842 


-3.62 


N821 


24.7±2.5 


2.47 


1.61 


196±26 


2.48 


-2.81 


1.18 


-3.13 


N1399 


20.5±1.6 


5.42 


3.60 


312±41 


59.7 


-1.79 


10.60 


-2.53 


N1600 


68.5±6.6 


19.14 


16.90 


307±40 


159 


-2.05 


9.83 


-3.24 


N2300 


27.0±2.6 


3.30 


3.40 


269±35 


23.3 


-2.16 


5.27 


-2.81 


N2778 


23.3±3.4 


0.557 


0.359 


171±22 






0.621 


-2.76 


N2832 


96.8±9.4 


14.9 


10.56 


349±45 


123 


-1.94 


18.00 


-2.77 


N3115 


9.8±0.6 


2.31 


1.59 


278±36 


4.74 


-2.53 


6.15 


-2.41 


N3377 


11.6±0.6 


0.891 


0.216 


131±17 


0.713 


-2.47 


0.176 


-3.09 


N3379 


10.8±0.7 


1.69 


0.822 


201±26 


4.29 


-2.29 


1.33 


-2.79 
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23.6±1.5 


2.51 


1.27 


206±27 


2.87 


-2.64 


1.49 


-2.93 
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31.7±6.2 


3.28 


2.22 


185±24 


10.4 


-2.36 


0.900 


-3.39 


N4278 


15.3±1.7 


1.90 


1.25 


270±35 


13.6 


-1.98 


5.36 


-2.37 


N4291 


26.9±4.1 


1.65 


1.11 
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17.5 


-1.81 


5.27 


-2.33 


N4467 


16.7±1.0 


0.0667 


0.0355 


87±11 






0.0256 


-3.14 


N4472 


16.7±1.0 


10.87 


8.94 


273±36 


28.5 


-2.51 


5.65 


-3.20 


N4473 


16.1±1.1 


1.85 


0.934 


188±25 


3.48 


-2.46 


0.971 


-2.98 


N4486 


16.7±1.0 


9.12 


9.04 


345±45 


38.6 


-2.36 


17.04 


-2.72 


N4486B 


16.7±1.0 


0.109 


0.0358 


178±23 
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-0.541 
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-1.68 


N4552 


15.7±1.2 


2.37 


1.56 


269±35 


4.79 
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5.27 


-2.47 


N4564 


14.9±1.2 


0.767 


0.416 


153±20 


2.48 


-2.24 


0.367 


-3.05 


N4594 


9.9±0.9 


5.10 


3.08 


248±32 


7.39 


-2.65 


3.59 


-2.93 


N4621 


18.6±1.9 


4.07 


2.34 


222±29 


3.39 


-2.84 


2.13 


-3.04 


N4636 


15. Oil. 1 


3.83 


3.16 


180±23 


2.22 


-3.15 


0.791 


-3.60 


N4649 


17.3±1.3 


7.85 


6.00 


331±43 


44.3 


-2.14 


14.02 


-2.63 


N4660 


13.2±1.3 


0.226 


0.118 


211±28 


2.41 


-1.70 


1.67 


-1.85 


N4874 


100.9±9.8 


26.1 


22.3 


230±30 


225 


-2.00 


2.51 


-3.95 


N4889 


91.6±8.9 


18.2 


11.9 


373±49 


265 


-1.68 


24.63 


-2.68 


N6166 


131±13 


28.2 


18.6 


311±41 


330 


-1.75 


10.45 


-3.25 


N7332 


23.5±2.3 


1.05 


0.193 


125±16 






0.141 


-3.14 


N7768 


114±11 


15.62 


9.79 


224±29 


101 


-2.03 


2.22 


-3.64 



NOTE.-Distances in Mpc. Lbuigo in 10"^Lq. Mbuigo in W^^Mq. CTc in km s"^ 
M, in IO^Mq. 
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Fig. 1. — Correlations between black hole mass and: (a) ^-band bulge luminosity; (b) bulge 
mass. Masses are in units of solar masses and luminosities in solar luminosities. Dashed lines are 

M,/ Mq = lO-^Lbuigc/i© (left panel) and M,/ M© = lO-^Mbuw/Mo (right panel). 
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Fig. 2. — Ratio of black hole mass computed by Magorrian et al. (1998), Mgt, to black hole mass 
computed from the M, — a relation, M,, as a function of the radius of influence of the nuclear BH. 



11 



o 



in 
d 



1 1 


1 


1 1 1 


1 1 1 1 


1 


- 




A'- 


\ 


- 


- 




/\i 


\ 


- 






/ \ 


\ 




- 




/ '\\ 


\ 
\ 


- 






/ ' w 


\ 




- 




/ ' \\ 


\ 


- 


- 




/ ' \ 


y \ 


- 


_ 




/ ' \ 


\ \ 


_ 






/ 1 


x\ \ 






\ 


1 -.-^ 1 1 


"y^^^-- ,- "",- 


__^^ 



-4 



-2 







log X 



Fig. 3. — Frequency function of logx where x = M,/Mt,uigc- Heavy solid line was derived from 
BH masses computed via the M, — a relation, equation (1); data are shown as the large dots. 
Dashed line was derived using the Magorrian et al. black hole masses; data are shown as the small 
dots. Thin solid line is the best-fit Gaussian approximation to A^(logx). Each curve assumes a 
measurement uncertainty in logx of 0.15. 
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Fig. 4. — (a) iV(logx) computed separately for the high-Mbuigc (thick hne) and low-Mbuigc (thin 
hne) galaxies in Table 1. (b) Effect of varying the assumed slope of the M, — a relation. Thick line: 
a = 5.25. Thin line: a = 3.75. Each curve assumes a measurement uncertainty in log a; of 0.15. 



